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Abstract 

We investigate the three-body Coulomb breakup of a two- neutron halo nucleus ^^Li. We use the 
coupled-channel ^Li + n + n three-body model, which includes the coupling between last neutron 
states and the various 2p-2h configurations in ^Li due to the tensor and pairing correlations. The 
three-body scattering states of ^^Li are described by using the combined methods of the complex 
scaling and the Lippmann-Schwinger equation. The calculated breakup cross section successfully 
reproduces the experiments. The large mixing of the s-state in the halo ground state of ^^Li is 
shown to play an important role in explanation of shape and strength of the breakup cross section. 
In addition, we predict the invariant mass spectra for binary subsystems of ^^Li. It is found that 
the two kinds of virtual s-states of ^Li-n and n-n systems in the final three-body states of ^-"^Li 
largely contribute to make low-lying peaks in the invariant mass spectra. On the other hand, in the 
present analysis, it is suggested that the contributions of the p-wave resonances of ^"^Li is hardly 
confirmed in the spectra. 

PACS numbers: 21.60.Gx, 21.10.Pc, 25.60.Gc, 27.20. +ii 
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I. INTRODUCTION 



Recent development of RI beam experiments reveals exotic and interesting properties 
of unstable nuclei beyond the stability line In particular, the neutron halo structure 
is one of the most interesting topics, and many studies have been performed from both 



experimental jl| and theoretical (2| sides since the discovery of the halo structure in ^^Li |3|]. 
The halo nuclei, such as ^He, ^^Li, ^^Be, and ^^Be, have been observed to have large matter 
radii [l| and narrow momentum distributions of the last neutrons ^|. These facts reflect the 
ground-state structure in halo nuclei, in which the weakly-bound last neutrons are spread 
widely far from the core nucleus. 

In addition to the exotic structure in the ground states, breakup reactions of halo nuclei 
are expected to provide us with much information on the excited states, which are mostly 
located above the particle thresholds. The transition mechanism into the excited states is 
an important issue to be understood. Experimentally, the Coulomb breakup reactions have 
been performed to search for the electric dipole responses of halo nuclei jsj. The observed 
cross sections commonly show the low-lying enhancement above the breakup thresholds 



121 , and this enhancement has been considered to be related to the exotic halo structure. 



For one-neutron halo cases such as ^^Be [6|, the mechanism of the Coulomb breakup reaction 
has been discussed as the direct breakup from the halo ground state to the non-resonant 
continuum states. The low-lying enhancement can be interpreted as the reflection of the 
weakly binding halo structure. On the other hand, for two-neutron halo cases such as ^He 
and ^^Li, the breakup mechanism is complicated because the two-neutron halo nuclei are 
broken up to the core + n + n three-body scattering states due to the Borromean nature 
and these final scattering states contain various kinds of correlations, such as of the binary 



subsystems of core-ra and n-n 



m]i5\. 



It has been observed in ^^Li that the amount of the (lsi/2)^ component (45±10 %) of the 
halo neutrons is comparable to that of This fact indicates the breaking of the 

N = 8 magic number in the ground state of ^^Li. For ^°Li, the virtual s-state is suggested 
from a large negative value of the observed scattering length [3]. The experiments on 



Based 



^"^Li also support the existence of p-wave resonances at low excitation energies 
on those observed properties of ^"^Li and ^^Li, it is necessary to understand the breakup 
mechanism of ^^Li using the Coulomb response. 



Theoretically, the exotic properties of Li have been studied by using several kinds of 



approaches, such as the ^Li + n + n three-body models 



variational method 



13 



191 ] ■ no-core shell model |20|, 



fermionic molecular dynamics |2l|, antisymmetrized molecular dynamics 



221. stochastic 



microscopic cluster model (2J], and so on. For the Coulomb breakup, 
three-body models have been often used to investigate the mechanism of the breakup into 
the three-body final states. Most of the studies based on the three-body models assume 
to make the lsi/2 orbit degenerated with the Opi/2 one energetically by using the different 
^Li-n interactions for even and odd parity states 
reproduce the observed large s-wave mixing in ^^Li 



25 



261]. However, those models cannot 



26| . Furthermore, Esbensen et al. 27] 



have mentioned that the ^Li + n + n three-body model fails in explanation of the observed 
charge radius and the dipole strength in ^^Li consistently. In order to overcome these 
difficulties in ^^Li, a new approach beyond the ^Li + n + n three-body model has been 
desired. 
Myo et al. 



28 



30| have attempted to understand the exotic properties of ^^Li by taking 
into account the tensor and pairing correlations in ^Li using the tensor-optimized shell model 



(TOSM) 3l|, l32|. In their approach, the s-wave virtual states and the p-wave resonances in 



^°Li are reproduced simultaneously by taking into account the 2p-2h excitations of the ^Li 
core coming from the tensor and pairing correlations. The Pauli principle between the last 
neutron and the ^Li core plays an important role. This dynamical description has also been 
applied to the ^^Li system by using the coupled-channel ^Li + n + n model. In this model, 
the tensor correlation produces the specific 2p-2h excitation from the Osi/2 orbit to the Opi/2 



orbit 



311 ]. and the Pauli-blocking between the excited Opi/2 neutron in the ^Li core and the 



last halo neutrons dynamically pushes up the energy of the Opi/2 orbit in ^^Li 28|, |29|]. As 
a result, the energies of the p-wave and s-wave states of ^^Li get close to each other, and a 
large s-wave mixing in the ^^Li ground state is brought. Furthermore, this coupled- channel 
^Li + n + n model is shown to reproduce various physical quantities such as the matter and 
charge radii, and the dipole strength of ^^Li consistently. 

In this study, we investigate the Coulomb breakup reaction of ^^Li by using the coupled- 



channel ^Li 



n 



n three-body model proposed in Refs. 28|, l29|] . The aim of this work is 



to investigate the role of the dynamical mixing of (lsi/2)^-configuration in the ground state 
on the reaction mechanism of the three-body Coulomb breakup of ^^Li. The other aim is 
to get the knowledge of the binary correlations of ^Li-n and n-n in the breakup final states. 



For these purpose, it is necessary to describe the three-body scattering states with a correct 
boundary condition. This problem is solved by applying the complex-scaled solutions of the 
Lippmann-Schwinger equation (CSLS) Q, [la]- This method has been recently developed 
and successfully applied to the Coulomb breakup reactions of ^He and the a + d scattering 
using the a + p + n model for ^Li jssl - 



The Coulomb breakup cross sections have been analyzed as functions of subsystem en- 
ergies to understand the correlations in two-neutron halo nuclei 



13 



15 



26 



34 



35|. For 



two-neutron halo nucleus ^He, the breakup mechanism of ^He has been analyzed based on 
the a + n + n three-body model, in which the employed a-n interaction reproduces the ob- 
served phase shifts but is not enough strong to produce the virtual s-states. In our previous 
works of ^He Q, [isl, it was found that the characteristic shape of the breakup cross section 
is weakly influenced by the halo structure in the ground state. The magnitude and the peak 
position of the cross section can be explained with a dominant contribution from the final 
state interactions (FSI) which produce the ^He(3/2~) resonance and the n-n virtual state. 



In contrast to ^He, the Coulomb breakup cross sections of ^^Li have been shown to 
depend strongly on the structure, in particular, on the s-wave mixing in the halo ground 



state 



13 



13 



26| , which predict a small 



The simple ^Li + n + n model calculations 
s-wave mixing (about 20 %), reproduce the low-energy enhancement in the cross section. 
However, the energy position of the enhancement is higher than the observed one [3]. In this 



paper, it is shown that this serious problem in the ^Li 



n 



n model is solved satisfactory. 



Furthermore, the previous ^Li + n + n model suggests that the virtual s-state and the pi/2 
resonances of ^"Li are observed in the strength distribution together with the n-n correlation 
as functions of subsystem energies j26|. These problems are also discussed in the present 
coupled- channel ^Li + n + n analysis. 



This paper is organized as follows. In Sec. HIl we expain the coupled-channel ^Li + n + n 
three-body model and the formalism of the Coulomb breakup reaction using CSLS. In 
Sec. mil we show the results of the Coulomb breakup cross section of ^^Li and the in- 
variant mass spectra of its binary subsystems, and discuss the breakup mechanism of ^^Li. 
All results and discussion are summarized in Sec. HVl 
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II. FRAMEWORK 



A. Coupled-channel three-body model for ^^Li 

We give a brief explanation of the coupled- channel ^Li + n + n three-body model of ^^Li 



employed here, the detail of which is given in Refs. 28|, |29|. To solve many-body correlations 
not only of the last two neutrons but also of the coupling with degrees of freedom in the ^Li 
core, we start with the Schrodinger equation: 

H^j. = E^j., (1) 

where H and $j7r are the total Hamiltonian and the wave function with spin-parity J'" for 
^^Li, respectively. 

The wave function $j7r is described as follows: 



^j. = Y, ^cA (^Li) ® x% (nn)] , (2) 



where A is an antisymmetrizer, and $^(^Li) and ^'^(rm) are the wave functions for the ^Li 
core and the last two neutrons, respectively. The index c represents the quantum numbers 
of the ^Li core configurations. The channel amplitude Cc is determined by solving the 
coupling between ^Li and the relative wave functions of the last two neutrons using the 



28 



29| for the configurations of the ^Li 



coupled-channel equation Eq. ([T]). Using TOSM 
core, we employ two kinds of 2p-2h configurations: (nn) j7r=o+ for the pairing correlation and 
{pn)j7r^i+ for the tensor correlation, in addition to the Op-Oh one of (05)^(0^3/2)^(0^3/2)^. 
Thus our model is considered to be an extension of the usual ^Li + n + n model, because 
this coupled-channel model is equivalent to the simple ^Li + n + n model when the only 
Op-Oh configuration is taken for ^Li. The Pauli principle between the last two neutrons is 
taken into account by antisymmetrizing the relative wave function x'^(nn). For the ^Li-n 
part, we take into account the Pauli principle by projecting out the Pauli forbidden states 
occupied by the ^Li core from the relative wave function x^(nra). 

Here, the total Hamiltonian for the coupled-channel ^Li + n + n three-body system is 
given as 

3 2 

H = J2f,-t^ + Vn-n + 5^ Kore-n(r^) + M^Li), (3) 
i=l i=l 
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FIG. 1. (Color online) Coordinate sets in hybrid- VT model for Li. 

where Tj and are kinetic energies of each cluster and the center-of-mass motion of the 
^Li + n + n three-body system, respectively. For the n-n interaction V^_„, we employ the 
realistic Argonne v8' force [s^]. The relative coordinate represents that between the ^Li 
core and the i-th last neutron, and Kore-n. is the ^Li-n interaction. For Kore-n, we employ 
a potential folding an effective A^A^ interaction with ^Li core density 38[. We here use the 
modified Hasegawa-Nagata (MHN) potential 39|, |40| obtained in the G-matrix calculation as 
the effective A^A^ interaction. We introduce one parameter 5, which enhances the strength 
of the intermediate range of MHN potential from the original one in order to describe 
the starting energy dependence dominantly coming from the tensor force in the G-matrix 



calculation [28|. In the present calculation, we take this parameter 6 as 0.1893 to reproduce 
the observed two-neutron separation energy of '^^Li. 

The microscopic internal Hamiltonian for the ^Li core, h{^Li), is given as 

9 9 

Here, ti and tcm are the kinetic energies for each nucleon and the center-of-mass motion 
in the ^Li core, respectively. The two-body A^A^ interaction Vij, whose detail is given in 
Ref. 3l|, consists of central, spin-orbit, tensor, and Coulomb terms. Using the Hamiltonian 
in Eq. @, the 2p-2h configurations in the ^Li core are involved by the pairing and tensor 
correlations. 

The wave function ^'^(rira) of the last two neutrons are solved by using a few-body tech- 
nique. Here we employ the variational approach called the hybrid- V^T model 411 • In this 
model, we expand the relative wave function x'^(nn) with two kinds of basis functions, shown 
in FIG. [U One is the cluster-orbital shell model (COSM; V-type) and the other is the ex- 



tended cluster model (ECM; T-type). These two kinds of basis functions are important 
for expressing the ^Li-n and n-n correlations simultaneously. Radial components of each 
relative wave function are expanded using the Gaussian basis functions 42| . The number of 
the basis functions is determined to reach the convergence of the numerical results. In the 
present calculation, we employ 18 basis functions for each coordinates, and the range of the 
Gaussian is taken up to 60 fm. 



B. Complex-scaled solutions of the Lippmann-Schwinger equation 

We describe the ^Li + n + n three-body scattering states of ^^Li by using CSLS jl^ to 
calculate the Coulomb breakup cross section. Before going into the formalism of CSLS, we 
briefly explain GSM 43|. In GSM, the relative coordinates shown in FIG. [1] are commonly 
transformed as follows: 

u{e)^u-\e) = ^e^', (5) 

where U {6) is a complex scaling operator with a scaling angle 6 being a real number. The 
coordinate $, represents the set of the two relative coordinates in the three-body system 
as used in FIG. [TJ Applying this transformation to the Hamiltonian H, we obtain the 
complex-scaled Schrodinger equation is given as 

where and <l>^ are the complex-scaled Hamiltonian and the complex-scaled wave function 
given as 

= U{e)HU^\e) (7) 

and 

$^ = U{e)^{^) = e(5*^)-^$(^e'^), (8) 

respectively. The factor, e'-^*^) -^^ comes from the Jacobian in the volume integral, and / = 2 
for the three-body systems. By solving the complex-scaled Schrodinger equation given in 
Eq. with a finite number of basis functions such as Gaussian, we obtain the eigenstates 
and energy eigenvalues of as {$^} and {E^} with a state index n, respectively. 

All the energy eigenvalues {E^} are obtained on a complex energy plane, governed by 



the ABG theorem 



44 



45(1, and their imaginary parts represent the outgoing boundary con- 



ditions. In FIG. [21 we show a schematic distribution of energy eigenvalues of the three-body 



lm(E) 




FIG. 2. (Color online) Schematic picture of energy eigenvalues of the three-body Borromean system 
in CSM. 



Borromean system in CSM. In CSM, the resonances of a many-body system are obtained as 
the isolated poles with the basis functions. On the other hand, the energy eigenvalues of 
continuum states are obtained on the 2^-rotated branch cuts starting from different thresh- 
olds of two- and three-body decay channels, such as ^°Li + n and ^Li + n + in the case 
of ^^Li. This classification of continuum states in CSM imposes that the outgoing boundary 
condition for each open channel is taken into account automatically by the imaginary parts 
of energy eigenvalues. Using the classification of continuum states in CSM, we can describe 
three-body scattering states without any explicit enforcement of boundary conditions. 

Furthermore, the complex-scaled eigenstates satisfy the extended completeness relation 
(ECR) Q,0, 

consisting of bound states, resonances, and rotated continua for the state 



index n as 



\ nl \ n\ 



(9) 



where {$^, form a set of biorthogonal states 46] . This relation is used when we describe 
the scattering states with the Lippmann-Schwinger equation. 

In CSLS, we start with the formal solution of the Lippmann-Schwinger equation given as 

1 



^(^)(k,K) = 0o(k,K) +lim 



^0o(k,K) 



(10) 



where (/'o(k, K) is the solution of the asymptotic Hamiltonian Hq with two relative momenta 
k and K in Jacobi coordinates of the three-body system. The total Hamiltonian H is the 
same as given in Eq. ([3]), and the interaction V is defined by subtracting the asymptotic 



Hamiltonian Hq from H . 

For the three-body breakup reaction of ^^Li, the asymptotic Hamiltonian Hq consists of 
the kinetic part of the ^Li + n + n system and the internal Hamiltonian of ^Li since no 
binary subsystems have bound states because of the Borromean condition for ^^Li. Hence, 
we can define Hq and its solution 0o as 

3 

Ho = h{^U) + Y,^,-f,^, (11) 

i=l 

ff„^„.(^«„ + _ + _j^„, (12) 

(r. R|,;.„(k, K)) = -i-*g.(9Li) ® e^ '+'K-f , (13) 

where /i and M are the reduced masses corresponding to k and K, respectively. The coor- 
dinates r and R are conjugate to the momenta k and K, respectively. Here, we assume that 
the ^Li core is in the ground state with energy of e^i^^ after the breakup. The interaction V 
in Eq. ( fTOj) is given as 

2 

^ = V;ore-n(ri) + Vn-n- (14) 

i=l 

We consider the incoming scattering states in bra-representation, which are used as the 
final states in the Coulomb breakup cross section of ^^Li. Assuming the Hermiticities of H 
and V , Eq. (fTOl) is rewritten as 

(vI/(-)(k,K)| = (0o(k,K)| + (0o(k,K)|\>lim \ . . (15) 

e^o E - H + ie 

In CSLS, we express the Green's function in Eq. ( ITSl) in terms of the complex-scaled Green's 
function. The complex-scaled Green's function with the outgoing boundary condition, 
Q^{E), is connected to the non-scaled Green's function Q{E) as follows: 

lim ^ I . = g{E) = u~\e)g\E)u{e). (le) 

E -H + ie 
The explicit form of Q^{E) is defined as 

where ECR defined in Eq. ([9]) is inserted. Using Eqs. f lT6|) and f lT7|) . we obtain the incoming 
scattering state \E'*^~-' in CSLS as 

(vI/(")(k,K)| = (0o(k,K)| 

+ 5f(0o(k,K)|Vt/-i(^)|<)-^-l^(<|.^|t/(^). ^^^^ 



n 



10 



It is noted that the scattering states in Eq. (fT8|) consist of two terms: The first term 
describes the non-interacting three-body continuum state, which is the same as given in 
Eq. f llOp . The second term contains all information of FSI. Using Eq. ( llSp . we can extract 
the effect of each component from the breakup cross section. 



C. Coulomb breakup cross section 

The Coulomb breakup reaction is considered to be dominated by the El transition. In 
the present calculation, we obtain the Coulomb breakup cross section using the El transition 



48 



49|. 



strength and the virtual photon number from the equivalent photon method 

In CSLS, we calculate the momentum distribution of the El transition strength from the 
^^Li ground state into the ^Li + n + n three-body scattering states. The distribution is 
given as 



(fB{El) 1 



dkdK 2 Jgs + 1 



[¥~\\.,K)\\0{Eimy . (19) 



where $gs and Jgs are the wave function and the total spin for the initial ground state of 
^^Li, respectively. The wave function \i/'^~)(k, K) is that for the three-body scattering state 
of ^Li + n + n with relative momenta, k and K, and here, is described by using CSLS as 
shown in Eq. (|T8i) . Using the coordinate sets shown in FIG. [H, the El operator 0{E1) is 
given as 

dm{El) = ^ei?Fi„(R) = ^ (riFi^(f + r^Y^mih)) , (20) 

where Yi^ is the spherical harmonics and m is the z component of the operator. 
Using Eq. f lT9|) . the two-dimensional energy distribution is defined as 



d€ide2 J J dkdK. 



2 1,2 \ / fc2 T^2 \ 

(21) 



X 6 \ ei — ] S {62 



2fi J \ 2M 

where ei and 62 are the relative energies of each subsystem, such as ^Li-n and n-n. Similarly, 
the total energy distribution of the El strength is given as 



dE J J dkdK. , , 

11 



TABLE I. Resonance energies and decay widths of resonances of ^''Li, measured from the 
^Li + n threshold. The scattering lengths are for the 5^/2 virtual states of ^^Li. The experimental 
data for a, is also listed. 



1+ 



1" 



Exp. 



Energy (keV) 
Width (keV) 
as (fm) 



275 
150 



506 
388 



-6.8 



-45.0 



Reference 
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Using Eqs. f l2T]) and fl22l) and the equivalent photon method, we can calculate the Coulomb 
breakup cross sections. The two-dimensional energy distribution and total energy distribu- 
tion of the cross sections are given as 



(Pa IQn 



deide2 9hc 



and 



da 16 



71 



deide2 
dB{El) 



(23) 



(24) 



dE dhc"^^' ^" dE ' 
respectively. Here, Nei{E^) is the virtual photon number with the photon energy E^, and 

E^ is given as 

E-^ = ei + 62 + S2n = E — £:9Li + 82^ (25) 

where S2n is the two-neutron separation energy of ^^Li. 

From Eq. fl23|) . the invariant mass spectra for binary subsystems such as of ^Li-n and n-n 
are given as 

da r ^ d'^a 
dei J deide2 
where £1 is the relative energy of the binary subsystem. 



III. RESULTS 



A. Properties of ^°Li and ^^Li 



111 



We here discuss the properties of ^°Li and ^^Li obtained in the present coupled-channel 
^Li + n + n model. In TABEL HI we show the ^°Li results. We obtain the p-wave doublet 



12 



TABLE II. Ground-state properties of ^^Li. See the text for details. 





TliGor. 


Exp. 




O 1 1 


O 1 0_l_tJ 




o.oy 


Q 1 o-un 1 fi^ 






3.53±0.06^ 






3.71±0.2q^ 


Rr-u ( fm) 


2.43 


2.467±0.037f. 






2.423±0.034£ 




44.0 


45±iq£ 


P{{Vz/2?){%) 


2.5 




P{{Vl/2?){%) 


46.9 




P{{d,/2?){%) 


3.1 




P{{d,/2?){%) 


1.7 





^ Reference 
^ Reference 
^ Reference 
Reference 
Reference 
^ Reference 
s Reference 



521 
531 
521 
521 
M\ 
551 
161 



resonances, 1"'" and 2+, constructed from the couphng between the spins of the last neutron 
and ^Li(3/2~). Two resonances are located at 275 keV and 506 keV measured from the 
^Li + n threshold with the decay widths of 150 keV and 388 keV for 1+ and 2"*", respectively. 
There are also the s-wave doublet states of 1~ and 2~ in ^°Li in the same reason of the spin 
coupling as the p-wave case. The scattering length for the 2~ state is obtained as —45.0 
fm, quite a large negative value, and shows a good agreement with the experimental one 
within the error bounds [l3|. The present large negative value of indicates the existence 
of the virtual s-states in ^'^Li(2~). It is noted that the previous value of is about —18 
fm in Ref. {29 1, which is negatively smaller than the present one, because the smaller 2n 
separation energy of ^^Li, suggested from the old experiment, was adopted to determine the 
^Li-n interaction. 



13 



We show the ground-state properties of ^^Li in TABLE [TTl The calculated matter and 
charge radii, i?m and -Rch, are obtained as 3.39 fm and 2.43 fm, respectively, and are in a 
good agreement with the observed data. We also calculate the probabilities of the partial 
wave components P((/j)^) of the last two neutrons. The probability of is very large 

and consistent with the observed data 16|, which is almost equal to that of the (pi/2)^ 



component. The large s-wave mixing generates the neutron halo structure and also indicates 
the breaking of the = 8 magic number. From the results, the present coupled- channel 
^Li + n + n model using TOSM for the ^Li core well reproduces the structures of ^°Li and 
^^Li, in particular, those related to the s-wave properties. Other ground-state properties, 
such as the quadrupole moment and the spin dipole moment of ^Li and ^^Li, have been 
discussed, and it was shown that the present coupled-channel ^Li + n + n three-body model 
also explains those observed quantities 0, [3^ . The detailed analysis from the experimental 



side has been summarized in Ref. 



2l|. 



In the previous works 28|, l29|] , it has been discussed that the 2p-2h configurations induced 
by the tensor and pairing correlations in the ^Li core are essential to explain the various 
physical observables in ^°Li and ^^Li. It was shown that in ^Li, the specific excitations of 
neutrons into the pi/2 orbit are occurred due to the tensor and pairing correlations, so that in 
^°Li and ^^Li, the coupling between the excited pi/2 neutron in ^Li and the last pi/2 neutrons 
give rise to the Pauli-blocking, which plays to lose the total energy of ^°Li and ^^Li with 
p-wave configurations. As a result, the s-wave configurations dynamically gain the energy 
and are largely mixed, which are sufficient to explain the s-wave properties of the two nuclei. 



It was found that t 



le tensor correlation gives the stronger effect of the Pauli-blocking than 



the pairing case (28 



29[. The large matter and charge radii of ^^Li can be the consequences 



of the Pauli-blocking effect and becomes the important quantities for the Coulomb breakup 
reaction, because the charge radius affects the El transition from the viewpoint of the sum 
rule value and this transition can be dominant in the Coulomb breakup reaction. Actually 



in Ref. 36| the low-energy El strength of ^^Li is sensitive to the s-wave component of the 
ground state. Esbensen et al. discussed j^^] the relation between the Coulomb breakup 
strength and the charge radius of ^^Li. They concluded that the simple ^Li + n + n model 
cannot explain the both quantities simultaneously. This problem is solved in the present 
coupled- channel ^Li + n + n model including the tensor and pairing correlations in the ^Li 
core. 



14 



To see the effect of the channel-couphng, we also performed the single channel calcula- 



tion [29| using the inert Li core with the Op-Oh configuration and adjusting Kore-n to fit the 
experimental value. The result shows that the (pi/2)^ configuration of last two neutrons 
dominates the ^^Li ground state with the probability as 90.6%, which preserves the p-shell 
magic number of neutrons. Instead, the configuration is mixed very small by only 

4.3%. Accordingly, the matter and charge radii are obtained as 2.99 fm and 2.34 fm, re- 



spectively 



291], both of which are smaller than the experimental values shown in TABLE HTl 



These results indicate that the halo structure is not so developed and the inert core assump- 
tion cannot explain the properties of ^^Li. For ^°Li, the scattering lengths of s-wave states 
do not show the negative value, which does not support the existence of the virtual s-state. 
The results of ^"^Li and ^^Li in the single channel calculations are significantly different from 
the results of the coupled ^Li + n + n model. This fact means that the single channel 
calculation using an inert core is inadequate to understand the structures of ^"^Li and ^^Li 
consistently. 



B. Coulomb breakup cross section 

We calculate the Coulomb breakup cross section of ^^Li using Eg. with the El 
strength distribution and taking care of the experimental resolution 5^. The target is Pb 
and the incident energy of the ^^Li projectile is 70 MeV/nucleon. The cross section measured 
from the ^Li + n + n threshold energy is shown in FIG. [31 It is found that the results show 
a good agreement with the experiment [l^ for shape and magnitude over the whole energy 
region. The distribution shows a low-lying enhancement at around 0.25 MeV and rapidly 
decreases as the energy increases. In the CSM calculation, there is no three-body dipole 
resonance, whose decay width can make a visible structure on the cross section. 

One of the characteristics of ^^Li is quite a large mixing as 45 % in the ground 

state which generates the halo structure. It is interesting to see the effect of the large s-wave 
mixing on the Coulomb breakup strength. We here calculate the El strength distributions 
using the different ^^Li wave functions, in which the coupling of the 2p-2h configuration 
involved only by the neutron pairing correlation in the ^Li core is taken into account. This 
restriction of the correlation in ^Li leads to the small component as 21.0 % in the 

ground state |36|. The calculated distribution is shown in FIG. [Has red (dashed) line. The 
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distribution shows a relatively small strength at the peak energy, the magnitude of which 
is about a half of the original one with a large s-wave mixing. The result indicates that a 
large s-wave mixing in the ground state plays a significant role in reproducing the observed 
low-lying enhancement in the breakup strength. This point has also been discussed in the 
previous analysis by changing the s-wave properties of ^^Li |36j. The explicit effect of the 
virtual s-state of ^°Li in the Coulomb breakup would be discussed using the invariant mass 
spectra in the next subsection. It is suggested that the large s-wave mixing enhances the 
dineutron correlation in the ^^Li ground state 29| . 

It is important to clarify the effect of the correlations in the ^Li core on the El strength 
distribution. For this purpose, we compare our coupled- channel calculation including only 
the pairing correlation, and that of the simple ^Li + n + n model assuming an inert ^Li 



core [^J, which gives 20.6 % of the component and —5.6 fm scattering length of 

the s-wave state of ^Li-n. The both wave functions contain almost the same amount of 
the s-wave component in the ^^Li ground state. In two kinds of results, the El strength 
distributions commonly have peaks at around 0.5 MeV, however, there exists the large 
difference of the strength at around the peak energy. This is due to the fact that about 15 % 
of the integrated strength in our calculation escapes to the higher excited ^^Li states having 
the excited components of the ^Li core. 




Energy [MeV] 

FIG. 3. (Color online) Coulomb breakup cross section of ^^Li, measured from the ^Li + n + n 
three-body breakup threshold. The red (solid) line represents the cross sections calculated with 
CSLS. The experimental data are taken from Re. jl^ . shown as the open squares with errorbars. 
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FIG. 4. (Color online) Comparison between the El strength distribution with different (si/2)^ 
components. The black (solid) line represents the present result used in FIG. [3l The red (dashed) 
line is the result using the ground-state wave function with (si/2)^ = 21 %. The blue (dotted) line 



is the result taken from Ref. 



26| 



From these comparisons, it is summarized that the large s-wave mixing in the initial 
ground state of ^^Li and the correlations in the ^Li core play the essential roles in reproducing 
the Coulomb breakup cross section, in particular, the position and the magnitude of the low- 
lying enhancement simultaneously. Furthermore, the large g-wave mixing gives the strong 
influence on the breakup mechanism. In our previous work 



we have discussed that the 



s-wave mixing in the ground state enhances the direct breakup process into the ^Li + n + n 
states instead of the sequential one via the p-wave resonances in ^''Li. In fact, when the 
s-wave mixing is about 40% in the ground state, the three-body direct breakup process 
exhausts the 66% of the integrated El strengths. This result is much different from the ®He 
case, in which the sequential process via the ^He(3/2^) resonance dominates the breakup 



reaction 



46| . This difference between ^^Li and ^He can be understood as the effect of two- 



neutron s-wave component on the Coulomb breakup process. It should be here noticed that 



the discussion on the breakup mechanism in Ref. |36| are based on the strength distribution 
calculated in CSM, in which the s-wave virtual states of ^°Li cannot be separated from the 
continuum states. Such a calculation cannot distinguish the sequential process via the s- 
wave virtual states from the direct breakup process. It is important to derive the invariant 
mass spectra of the ^Li + n subsystem in order to estimate appropriately the contribution 
of the sequential breakup process via the virtual state in the final states, which would be 
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FIG. 5. (Color online) Comparison between the calculated Coulomb breakup cross sections of ^^Li. 
The black (solid) line represents the same result in FIG. [3l The red (dashed) line is the result 
without FSI. 



shown in the next subsection. 

In addition to the ground-state properties of ^^Li, it is also interesting to see the effect of 
FSI on the Coulomb breakup cross section. The FSI is defined in Eq. ( fT4l) . In CSLS, all the 
effects of FSI in the scattering wave functions are included in the second term in Eq. 0181) . 
We can drop off the second term in the calculation of the cross section to examine the effect 
of FSI, while we do not change the initial ground-state wave function. The cross section 
without FSI is shown in FIG. [5] as red (dashed) line and has a broad peak structure at 0.5 
MeV. It is found that the magnitude is much smaller than the full results including FSI, while 
there is a very small difference between two results in the higher energy region above 1 MeV. 
From this analysis, it is concluded that in addition to the initial state properties of ^^Li, FSI 
gives a significant effect to create the low-lying enhancement in the Coulomb breakup cross 
section of ^^Li. This large effect of FSI has been shown in the previous analyses 
for the Coulomb breakups of ^^Li and ^He. 
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C. Invariant mass spectra of binary subsystems of Li 

To see the effect of FSI in more detail, we calculate the invariant mass spectra using 
Eq. ( 126|) in the Coulomb breakup reaction of ^^Li. In FIG.[6l we show the results as functions 
of the relative energies of ^Li-n and n-n subsystems as panels (a) and (b), respectively, 
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together with the results calculated without FSI. It is found that both spectra have sharp 
peak structures commonly below 0.1 MeV. From those results, the peaks in the invariant 
mass spectra are understood to come from FSI. In FIG. M (b), the peak seen in the n- 
n invariant mass spectra is caused obviously by the n-n virtual state, because FSI, Vn-n, 
produces no resonance. Such a peak due to the n-n virtual state is also seen in the Coulomb 
breakup reaction of ^He 15 1. 



For the invariant mass spectrum of the ^Li-n subsystem, Hagino et al. 26[ discussed that 
both the s-wave virtual state and the p-wave resonance of ^"Li contribute to the El strength 



distribution using the simple ^Li 



n 



n three-body model without the tensor correlation 



in the ^Li core. On the other hand, as shown in FIG. [6] (a), our prediction indicates only 
a single prominent peak below 0.1 MeV. This peak is considered to come from the s-wave 
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FIG. 6. (Color online) Invariant mass spectra for ^Li-n and n-n binary subsystems. The panels 
(a) and (b) represent the results for ^Li-n and n-n subsystems, respectively. The red (solid) lines 
show the results with FSI and the black (dashed) ones are those without FSI. 
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FIG. 7. (Color online) s-wave and p-wave components of the invariant mass spectra for ^Li-n. 
The red (solid) and blue (dashed) lines represent the results for relative s- and p-wave components, 
respectively. The black (thin) line is same as the result with FSI in FIG. [6] (a). Two arrows indicate 
the positions of resonance energies of 1"*^ and 2"*". 



virtual state due to Kore-n- 

To see explicitly the partial-wave contributions in the ^Li-n invariant mass spectrum, we 
calculate the decomposed spectra as shown in FIG. [71 The results show that the s-wave 
component has a peak below 0.1 MeV, which comes from the virtual s-state of ^°Li. The 
p-wave component has a broad bump at around 0.15 MeV. It is found that the bump energy 
does not correspond to the p-wave resonance energies in ^"Li indicated by two arrows in 
FIG. [7l This result means that the p-wave contribution in the spectra comes from the non- 
resonant continuum states of ^Li-n. Thus we can conclude that the shape of the invariant 
mass spectra of ^Li-n is mainly determined by the virtual s-state in ^°Li, while the non- 
resonant p-wave contribution gives some amount in the spectra, which becomes dominant at 



higher energies than 0.2 MeV. This conclusion contradicts the result in Ref. |26|, in which 
the p-wave resonance has a sizable contribution to the strength and the virtual s-state seems 
to give a small effect on the strength in comparison with the present result. 

The reason why the p-wave resonances are not observed in the present ^Li-n invariant 
mass spectra can be understood as follows; The p-wave resonances of ^°Li are located at 
0.275 MeV and 0.506 MeV for l"*" and 2"*", respectively, as shown in TABLE [B On the other 
hand, the breakup cross section has a peak at around 0.25 MeV as shown in FIG. [31 This 
peak energy is lower than the energies of the p-wave resonances of ^"^Li. The relation of 



20 



energies implies that the sequential breakup process via the p-wave resonances of ^"^Li is 
energetically not favored at around the peak energy of the cross section. From this energy 
condition, the p-wave resonances give a minor contribution to the Coulomb breakup reaction 



of ^^Li. In fact, since the observed breakup cross section shows a peak at 0.25 MeV 12|, it 
seems to be difficult to observe the p-wave resonances in the ^Li-n invariant mass spectra 
if those energies are higher than 0.25 MeV. On the other hand, the breakup cross section 



calculated by Hagino et al. 



26| has a peak at around 0.5 MeV, locating higher than the 



p-wave resonance energies of ^°Li. In that case, the sequential breakup via the p-wave 
resonances of ^°Li is favorably allowed, and can make the peak in the strength. 

The contributions of the s- and p-wave components in the ^Li-n invariant mass spectra 
can be seen more clearly in the two-dimensional energy distributions of the cross section, 
d?a / deide2i shown in FIG. El Here, Si and 62 are the relative energy of the ^Li-n subsystem 
and the energy of the relative motion between the center-of-mass of the ^Li-n subsystem 
and the other neutron, respectively. The distibutions are calculated by using Eq. f l2^ . It 
is found that the total distribution shown in FIG. [8] (a) has two peaks at small energies of 
El and 62-, respectively. The origin of those structures are confirmed by decomposing the 
distribution into the s- and p-wave components for the ^Li-ra states. In FIG. [8] (b), the 
s-wave component is shown and the concentration of the strength is found at ei = 0.03 MeV 
and £2 = 0.13 MeV. This peak near the zero-energy of ei is a consequence of the virtual 
s-state of the ^Li-n system. In FIG. [8] (c), the p-wave component shows a peak at ei = 0.12 
MeV and 62 = 0.03 MeV and the distribution is relatively wider for ei than the s-wave 
case, which is consistent to the non-resonant breakup of the p-wave component in the ^Li-n 
system. It is also found that the s- and p-wave components have the similar magnitudes. 
Small bumps are obtained at small energies of £2 in FIG. El (b) and of ei in FIG. El (c) for 
the s- and p-wave components, respectively. This result suggests a coupling between s- and 
p-wave components. 



IV. SUMMARY 

We have investigated the Coulomb breakup reaction of two-neutron halo nucleus ^^Li 
using the coupled-channel ^Li + n + n three-body model. The three-body scattering states 
of ^^Li are described in terms of the complex-scaled solutions of the Lippmann-Schwinger 
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FIG. 8. (Color online) Two-dimensional energy distributions of the Coulomb breakup cross section 
of ^^Li, where ei and 82 are the relative energy of the ^Li-n subsystem and the energy of the relative 
motion between the center-of-mass of the ^Li-n subsystem and the other neutron, respectively. The 
panel (a) presents the total distribution, and the panels (b) and (c) present the results for s- and 
p-wave components for the relative motion of ^Li-n, respectively. 



equation (CSLS). In the present model of ^^Li, we take into account 2p-2h configurations 
of the ^Li core, which describes the tensor and pairing correlations in the ^Li core on the 
basis of the tensor-optimized shell model. As a result, we can reproduce the breaking of 
the = 8 magic number due to the large s-wave mixing, which brings the neutron halo 
structure in ^^Li, and also explains the s-wave property of ^°Li simultaneously. In this paper, 
we calculated the transition from the ^^Li ground state into the three-body scattering states 
using CSLS by the Coulomb response. In CSLS, the Green's function using the complex- 
scaled wave functions provides with the three-body scattering states of ^^Li, which include 
much information of the correlations not only of three-body system but also of the binary 
subsystems of ^Li-n and n-n. 

The calculated Coulomb breakup cross section of ^^Li into the three-body scattering states 
shows the low-lying enhancement and well reproduces the observed data over whole energy 
region. The magnitude of the low-lying enhancement is sensitive to the s-wave mixing in 
the initial ^^Li ground state, the value of which is obtained as 44% in the present calculation 
and reproduces the experimental value. The excitation of the ^Li core is also important 
to reproduce the total breakup strength, because about 15% of the strength escapes to the 
higher energy region as the component of the core excitation in the present coupled- channel 
approach. This component reduces the low energy strength corresponding to the experiment. 
These two effects nicely work to determine the distribution of the cross section of ^^Li. 

It is confirmed that the low-lying enhancement in the cross section is also affected by 
the final-state interactions (FSI) in the ^Li + n + n three-body system. This result is a 
similar feature to the ^He breakup case except for the large s-wave mixing in the ground 
state. To determine what kinds of FSI dominate the Coulomb breakup reaction of ^^Li, we 
calculate the invariant mass spectra as functions of the energies of the two kinds of the binary 
subsystems of ^°Li as ^Li-n, and n-n. The calculated spectra show the sharp peak structures 
below 0.1 MeV from the thresholds in both cases. It is found that those peaks come from 
the virtual-state correlations in ^°Li and n-n subsystems in the final states, respectively. 
On the other hand, it is found that the p-wave resonances of ^"^Li give a minor effect on 
the breakup cross section. This is because the p-wave resonances of ^°Li are energetically 
located higher than the energy of the low-lying enhancement in the breakup cross section, so 
that the sequential breakup process via the p-wave resonances is mostly forbidden in energy. 
In relation to this fact, the effect of FSI on the p-wave component of the ^Li-n invariant 
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mass spectra is suggested to be small. This feature of the ^^Li breakup is very different from 
the ^He case, which shows a large effect of the p-wave component of the ^He-n states, and 
instead, a minor effect of the s-wave component. 
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